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ABSTRACT 

We investigate the transition from primordial Pop III star formation to normal Pop II star formation in the first 
galaxies using new cosmological hydrodynamic simulations. We find that while the first stars seed their host 
galaxies with metals, they cannot sustain significant outflows to enrich the intergalactic medium, even assuming 
a top-heavy initial mass function. This means that Pop III star formation could potentially continue until z ~ 6 
in different unenriched regions of the universe, before being ultimately shut off by cosmic reionization. Within 
an individual galaxy, Pop II stars overtake Pop III stars in 20-200 Myr, based on the amount of stellar feedback 
and metal production. 

Keywords: galaxies: formation — galaxies: evolution — methods: numerical — stars: formation — cosmol- 
ogy: theory 
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1. INTRODUCTION 

The first stars in the universe formed in gas devoid of met- 
als. This exotic environment may have caused the initial mass 
function (IMF) for Population III stars to be different from the 
modern day case. Namely, the high Jeans mass of metal- free 
gas suggests a top-heayv IMF (lAbel et al.ll2000t iBromm et al.l 
ll999tlAbel et al.ll2002tlYoshida et al.ll2003l) . In turn, the feed- 
back processes in the first stars may have been more drastic, 
promp ting the release of extrem e amounts of ionizing radi- 
ation ( Tumlinson & S hull 2000) a nd the occurrence o f pair 
instability supernovae (PISNe) (IHeger & Wooslevll2002l) . 

To explore the effects of these stars on their host galaxies, 
we developed a model for Pop III star formation and feed- 
back and implemented it into the adaptive refinement tree 
(ART) code, as described in a companion paper, Muratov et 
al. (2012a, hereby Paper I). Pop III stars were modeled to 
form in gas that was dense, partially molecular, and of primor- 
dial composition. Pop III SNe and ionizing radiation feedback 
were enhanced relative to their Pop II counterparts, and the 
first PISNe seeded the ISM with metals. We ran a suite of cos- 
mological simulations with this model, and found that the dy- 
namical impact of Pop III feedback depended strongly on the 
mass of t he galaxy. In agreement with previous work in the 
field (e .g. IBromm et al1l2003h IWhalen et al1l2008t IWise et all 
l2012ah . we found that PISNe were able to efficiently expel gas 
and metals fro m the Mh -lO^M^ h alos expected to host the 
very first stars (Teg mark et al.l[l997h . However, these effects 
were often temporary, as cosmological inflows of fresh gas 
restored them to the universal baryon fraction. The metals, 
which had previously escaped past the virial radius, also typi- 
cally fell back into the growing potential wells of the accreting 
galaxies, leaving the intergalactic medium (IGM) mostly pris- 
tine. In galaxies with mass Mi, > 10 7 M Q , most gas remained 
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bound even after a PISN event, and metals were not ejected 
past the virial radius. 

Since Pop III stars by definition only form in primordial 
gas, the large amount of metals released in P ISNe leads to 
the 's elf-termination' of Pop III star formation (Yoshi da et al.l 
2004). According to our findings in Paper I, this self- 
termination can only be local, as enrichment of the IGM and 
external halos is rather minimal from single PISNe. There- 
fore, determining the epoch when Pop III termination be- 
comes universal is a somewhat different question. Pop III star 
formation could be relevant for a much longer phase of cos- 
mic history if a Pop III star formed in every pristine halo with 
M~> 10 6 - 10 s M , as the abundance of such halos increases 
considerably with cosmic time. 

Population II star formation can commence in galaxies once 
they are either sufficiently massive to enable the rapid gas 
cooling by atomic hydrogen lines, or enriched enough to en- 
able efficient metal cooling ( Ostri ker & Gnedin] 1 1 9961) . Be- 
cause the feedback of Pop II stars is weaker than that of Pop 
III, Pop II star formation should ramp up rapidly in the host 
galaxy, provided that accretion from filaments continues to 
bring new supply. However, the relative weakness of the 
feedback, taken in conjunction with the plethora of primor- 
dial sites where Pop III stars may still form, means that the 
host galaxy, as well as the universe as a whole, are still in- 
fluenced by Pop III stars for some time after Pop II star for- 
mation begins. Understanding this transition quantitatively is 
relevant to whether future observational facilities such as the 
James Webb Space Telescope (JWST) will be able to observe 
galaxies in the Pop III dominance phase. Studies thus far have 
shown that the first galaxies generally sit on the brink of de- 
tectability by JWST dPawlik et alJl2mll20T2T) . 

In this paper, we follow the evolution of the galaxies de- 
scribed in Paper I through the epoch of dominance of the 
first stars. We study the transition from Pop III to Pop II 
star formation, and quantify the duration of this epoch. We 
also explore the effect of cosmic variance, and determine the 
prevalence and importance of Pop III stars at various cosmic 
epochs. 

2. SIMULATIONS 

A full description of our simulation setup, including the de- 
tails of both the Pop III and Pop II star formation recipes, 
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is presented in Paper I. Here, we outline the setup only 
briefly. We perform the simulations with the Eulerian 
gasdynamics+N-body adaptive refin ement tree (ART) code 
(iKravtsov et all 1997tlKravtsovll 9991 r2003HRudd et alJ2008t 
Gnedin & Kravtsovll201 ll) . We use a 256 3 initial grid with up 
to 8 additional levels of refinement. For most of our runs, we 
apply this grid to a 1 h Mpc comoving box and the WMAP- 
7 cosmology (Q m = 0.28, fl A = 0.72, h = 0.7, cr 8 = 0.817, 
VL), = 0.046, VIdm = 0.234). This gives us a DM particle mass 
of rn-DM = 5.53 x 10 3 M Q and a minimum cell size of 22 co- 
moving pc. We also employ a 0.5 ft -1 Mpc comoving box, 
where using the same grid, the DM particle mass is set to 
(Udm = 690 M and the minimal cell size is 1 1 comoving pc. 

Pop III stars formed in the almost pristine gas with the 
abundance of heavy elements below t he critical metallicity 
log 1() Z/Z =-3.5 (Bro mm et alJl20 01). Cells were allowed 
to form Pop III stars if the gas density exceeded a threshold 
nH,min, and a molecular hydrogen fraction threshold ///,.„„„. 
Through a series of convergence tests, we found that = 
10 4 cm" 3 and //j 2 , mi „ = 10~ 3 were appropriate values for the 
two thresholds. 

The Pop III prescription was designed to test the maximum 
possible effect of feedback, relying on an IMF that was top- 
heavy. Half of the Pop III stars formed as 170M Q parti- 
cles and were set to explode in PISNe. Each PISN injected 
27 x 10 erg of thermal energ y and 80 M of metals into the 
ISM dHeger & Wooslevll2002l) . The remaining 50% of Pop III 
stars formed as 1OOM particles that explode in type II SNe, 
generating 2 x 10 51 erg of energy. All Pop III stellar particles 
emited a factor of 10 more ionizing photons per second than 
their Pop II count erparts of the same mass (Schaerer 2002; 
I Wise & Cenll2009l) . A suite of cosmological simulations per- 
formed with this model revealed that Pop III stars drastically 
affected halos with M h < 3 x 10 6 M Q , but not halos of higher 
masses. Extended convergence tests revealed that without suf- 
ficient mass resolution, it was easy to miss these important 
dynamical effects. 

In gas that is enriched beyond the critical metallicity, Pop II 
star formation is prescribed according to the molecular-based 
star formation recipe that is pr esented in lGnedin et al.l (2009) 
and lGnedin & Kravtsovl (120111) . 

In Table Q] we list all of the simulations from the suite 
which we analyze in this paper. Some of the simulations 
that were used in Paper I for convergence tests and for de- 
termining the ideal parameters for Pop III star formation are 
not included in this study. Here, we primarily represent the 
three 1 ft -1 Mpc boxes with runs UnderDense-_nHle4_fid, 
UnderDense+_nHle4_fid, and OverDense_nHle4_fid. Ad- 
ditionally, we study runs with extreme feedback, where we 
increased PISNe energy (run OverDense_ExtremeSN) and 
ionizing photon emission (run OverDense_ExtremeRad) by 
a factor of 10 relative to the fiducial model. We also in- 
clude a run with a conservative "low-mass" Pop III IMF 
that assigns the same feedback parameters to Pop III as for 
Pop II stellar particles (run OverDense_LowMass) to ac- 
count for the present uncertainty in the Pop III IMF (e.g. 
lO'Shea & Normanll2007HGreif et al.ll20ll . 

In Paper I, we demonstrated that mass resolution was cru- 
cial for capturing Pop III star formation in halos of mass 
M/,~1O 6 M , and that halos with M h < 3 x 1O 6 M were most 
susceptible to Pop III feedback. Our fiducial 1 h~ l Mpc runs 
lacked the resolution to study these objects effectively, but 
using the same grid on a smaller box naturally increased the 
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Figure 1. Cumulative number of halos vs. mass at two epochs for the three 1 
h Mpc boxes and one H Mpc box. The smallest mass plotted corresponds to 
the earliest halo to form a Pop III star among all of our runs. Box OverDense 
(black) has an order of magnitude more halos than Box UnderDense+ (red) 
and Box UnderDense- (blue) across almost the entire range of masses con- 
sidered. While Box UnderDense+ and Box UnderDense- have similar DC 
mode values, the former clearly has more massive halos. One Box Under- 
Dense+ galaxy in particular is on par with the most massive galaxies in Box 
OverDense. The H Mpc box (green, ran OverDense_HMpc_HiRes) repre- 
sents only an eighth of the volume of the other boxes, and therefore hosts 
fewer massive galaxies while still representing an overdense region. 



resolution while keeping the computational cost of the simu- 
lations down. For this reason, we also use a 0.5 h~ l Mpc box 
(referred to as H Mpc) for run OverDense_HMpc_HiRes, to 
test the validity of our results in the low-mass regime. 

3. RESULTS 
3.1. Cosmic Variance 

Using the D C mode formalism for the gene ration of ini- 
tial conditions (Sirko 2005; Gnedi n~et alj|201 ll) allows us to 
test several representative regions of the universe without sac- 
rificing resolution, as would be needed were we to simulate 
a larger cosmological volume. With a single parameter that 
stays constant over time in a given simulation box, Ape, we 
encode the amplitude of density fluctuations on the funda- 
mental scale of the box. Many studies have been done to 
understand the effects of cosmic varian ce on the dark mat- 
ter halo mass function dTinker et al.N2.Q08h , however studying 
this effect convolved with hydrodynamics has thus far proven 
to be rather difficult. Here we present an analysis of the vari- 
ance of the three different 1 h~ l Mpc simulation boxes in our 
study. For reference, the DC mode values in the three boxes 
are A DC = -2.57,-3.35, and 4.04 in Box UnderDense-, Un- 
derDense+, and OverDense respectively. The H Mpc box has 
a DC mode of Adc = 5.04. Large, positive values of Aqc 
indicate an overdense region. 

The first Pop III stars form around scale factor a « 0.047 
(z ~ 20.3) in both Box UnderDense+ and Box OverDense. On 
the other hand, Box UnderDense- stalled significantly, and no 
star formation occurs until a m 0.073 (z ~ 12.7), translating to 
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Table 1 

Simulation runs 



Run 


Base grid 




dx (pc) 


m DM (M Q ) 


n H .mm(cm 3 ) 


Description 


Fiducial runs 


UnderDense— _nH 1 e4_fid 


256 3 


8 


22 


5500 


10000 


Underdense box, no high-mass galaxies, fiducial parameters 


UnderDense+_nHle4_fid 


256 3 


8 


22 


5500 


10000 


Underdense box, one high-mass galaxy, fiducial parameters 


OverDense_nH 1 e4_fid 


256 3 


8 


22 


5500 


10000 


Overdense box, fiducial parameters 


Increased mass resolution 














OverDense_HiRes 


512 3 


7 


22 


690 


10000 


High mass resolution 


OverDense_HMpc_HiRes 


256 3 


8 


11 


690 


10000 


0.5 /T 1 Mpc box, high spatial resolution 


Alternative feedback 














OverDense_ExtremeSN 


256 3 


8 


22 


5500 


10000 


Extreme PISNe (Sectionl3!4l 


OverDense_ExtremeRad 


256 3 


8 


22 


5500 


10000 


Extreme Pop III radiation field (Section[3.4) 


OverDense_LowMass 


256 3 


8 


22 


5500 


10000 


Pop III IMF and feedback mirror Pop II (Section[3.5) 



Column 1 .) Name of the run; 2.) Base grid, number of DM particles, number of root cells; 3.) Maximum number of additional levels of refinement; 4.) Minimum 
cell size at the highest level of refinement in comoving pc; 5.) DM particle mass in Mq; 6.) Minimum H number density for Pop III star formation in cm -3 ; 7.) 
Further description of the run. 



a time difference of 165 Myr. This wide range demonstrates 
immediately that the epoch when Pop III star formation be- 
gins is a strong function of local overdensity. If voids like 
the one represented by Box UnderDense- are indeed not en- 
riched by external sources, Pop III stars could have existed in 
these voids until the end o f cosmic reionization , at epochs that 
will be probed by JWST dHummel et al.ll2012l) and the La rge 
Synoptic Survey Telescope (LSST) (see lTrenti et al. 2009). 

Figure[T]shows the number of halos capable of hosting Pop 
III stars. Box OverDense clearly dominates over the other 
two across almost the entire mass range, and as halo mass 
strongly correlates with the density of central gas cells, Box 
OverDense is host to many more star forming galaxies. Box 
UnderDense- and Box UnderDense+ have similar numbers 
of low-mass halos with M < 10 7 M . However, the two sim- 
ulations very clearly diverge at the high-mass end. The most 
massive halo in Box UnderDense- has about 10 analogs in 
Box UnderDense+. The H Mpc box contains fewer halos 
per mass bin than Box UnderDense+, but per volume con- 
tains a higher density of massive halos, consistent with its 
higher Aoc value. At the time of formation of the first stars 
(around a = 0.05), none of the halos are more massive than 
2 x 1O 7 M . 

Figure [2] presents the star formation rate (SFR) density in 
each box for the runs with fiducial parameters. We can im- 
mediately see that SFRs vary by orders of magnitude among 
the three boxes. The most drastic aspect of this comparison is 
that Box UnderDense- does not form any stars until a = 0.073 
(z = 12.7). Pop II stars take a relatively long time to appear 
in the run, and the total mass of Pop II is only 2300 Mq by 
a = 0.1 (z = 9). Only two galaxies are able to form Pop III 
stars in this run. Such a narrow margin of error suggests that it 
is quite possible that with slightly different parameters for the 
initial overdensity, galaxies in this box may have never been 
able to form any Pop III stars before all halos were stripped 
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Figure 2. The SFR density vs scale factor for runs UnderDense-_nHle4_fid 
(blue), UnderDense+_nHle4_fid (red), OverDense_nHle4_fid (black), and 
OverDense_HMpc_HiRes (green). Global SFRs are significantly affected by 
the initial conditions in their respective simulation boxes. Pop III SFRs are 
denoted by open circles, Pop II SFRs by filled circles. 



of gas by external sources of ionizing radiation. 

Box UnderDense+ and Box OverDense have similar global 
SFRs at a = 0.065-0.08 (z = 14.4- 1 1 .5), as each box is being 
dominated by only a few galaxies at early times. However, the 
disparity between the mean overdensity of the boxes soon be- 
comes evident as Box OverDense becomes filled with a pop- 
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ulation of massive halos able to host dense cores in which 
Pop III, and subsequently, Pop II star formation is initiated. 
Box UnderDense+ continues hosting only one such galaxy 
for a large fraction of the duration of the simulation. Pop 
III star formation is completely stopped in Box UnderDense+ 
between a = 0.06 and a = 0.085, until an external halo of pri- 
mordial composition reaches the Pop III threshold density and 
later merges with the central galaxy. 

Using the many realizations of Box OverDense presented 
in Paper I, we find the variance in the SFR between realiza- 
tions within the same box to be low. Differences between 
the setup of our runs primarily affect Pop III stars, but at 
a = 0.055, the total number of Pop III stars did not depend 
on the gas density threshold we used for their formation. Fur- 
thermore, Pop II stars constitute most of the star formation 
past a = 0.07. Once a given galaxy transitions to Pop II as the 
dominant stellar population, none of the parametric variations 
for Pop III star formation affect the history of that galaxy. 
The Pop III SFR stays relatively constant between 10~ 5 and 
10~ 4 M Q yr" 1 throughout the entire simulation in Box Over- 
Dense. The range of values and relative constancy of this SFR 
up to z~6 are similar to the results found in l arger-scale (4 
Mpc) SPH simulations bv lJohnson et al.1 (120121) . 

The H Mpc box, represented here by run Over- 
Dense_HMpc_HiRes, samples a very overdense region of the 
universe. The Pop III SFR density in this run significantly ex- 
ceeds the others before a = 0.065. By a = 0.075, Pop II star 
formation has also become dominant in this box, but the to- 
tal SFR density has not reached the corresponding values for 
Box OverDense and Box UnderDense+. This disparity results 
from the fact that the 1 Mpc boxes ultimately sample more 
massive galaxies which are able to sustain high Pop II star 
formation rates. In addition, a larger fraction of galaxies sim- 
ulated in the H Mpc box have considerable gas blowout after 
PISNe, delaying Pop II star formation for a longer period of 
time. It follows that a halo from the H Mpc run typically has 
a lower SFR compared to one of equivalent mass from the 1 
Mpc runs at early times. 

3.2. Transition to Normal Star Formation 

To understand the environment of Pop III star formation in 
their host galaxies, we focus on run OverDense_nHle4_fid 
during the first 140 Myr after the first star forms. This cor- 
responds to the range of scale factor 0.048 < a < 0.07, or 
equivalently 20 > z > 13. After that epoch the most mas- 
sive galaxies in this box all form exclusively Pop II stars, and 
quickly build up a dominant mass of Pop II stars. Even more 
importantly, by a = 0.07, all Pop III stars in these galaxies 
have already exploded as SNe, and no longer provide radia- 
tive feedback. We w ill quantify the duration of the Pop III 
epoch in Section [3~3l 

Figure [3] shows the distance of each stellar particle from 
the center of its host galaxy at a = 0.07 vs the cosmic time at 
which the stellar particle was formed. Five of the most mas- 
sive star-forming galaxies are stacked together in this plot, 
showing the similarity in the assembly history of these galax- 
ies. We see that each of the five galaxies does not host more 
than a few Pop III stars. If a 100M Q star is the first to form in 
a given galaxy, it will not release any metals, and further Pop 
III star formation will be possible after the gas recollapses. If 
on the other hand a 170M Q star is formed, Pop III star forma- 
tion is rarely possible in the same galaxy again. The metals 
generated by a single PISN from a 170M Q star are enough to 
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Figure 3. Bottom panel: distance of each star from its host galaxy's cen- 
ter at a = 0.07 (z = 13) vs the cosmic time when the star formed, for five of 
the most massive star-forming galaxies run OverDense_nHle4_fid. Pop II 
stars are red points, the remnants of 100 Mq Pop III stars are blue symbols, 
and 170 Mq Pop III remnant tracer particles are black symbols. Individual 
galaxies are differentiated by different symbols. Every galaxy has a similar 
history: a few Pop III stars form and temporarily quench further star for- 
mation. However, eventually gas recollapses and self- sustained Pop II star 
formation begins within a 10-30 pc core, enriching the dense gas enough for 
Pop III star formation to be no longer possible within the galaxy. Pop III 
remnants often end up somewhat displaced from the Pop II core by as much 
as 100 pc. Top panel: the sum of SFRs for Pop III (solid) and Pop II (dotted) 
for these five galaxies, in black. For comparison, these rates are also plotted 
for run OverDense_LowMass, in red. 



switch the galaxy's primary mode of star formation to be Pop 
II by pushing the metallicity of gas above the critical threshold 
of l og ln Z/Z(7) = -3.5 (as has been previously demonstrated 
e.g. I Wise et al]|2012al) . A few exceptions to this pattern hap- 
pen when the sites of Pop III star formation within the galaxy 
are sufficiently spread apart (i.e. two Pop III stars form in 
two different halos that subsequently merge) or alternatively, 
the Pop III stars form sufficiently close together in time. The 
latter case is possible when many cells within the same giant 
molecular cloud reach the threshold density for star formation 
and several Pop III stars form before their feedback disperses 
the remaining gas in the cloud. This produces Pop III star 
multiples that are very close to each other in location and age. 
The spatial distributions of Pop III stellar remnants and Pop 

II stars within their host galaxy are significantly different. Pop 

III stellar remnants can be found out to 70 pc from the center, 
while most Pop II stellar particles at this epoch are within 20 
pc. Note that the centers of the galaxies are defined to only 
-5-15 pc, because of the limited number of DM particles and 
force resolution. More massive galaxies have better-defined 



The epoch of Population III stars 



5 



centers. Both stellar populations generally form close to the 
center of the host galaxy, but Pop III stars typically form at 
earlier times when the galaxy is not very massive, and the 
spatial concentration of DM particles is not very high. In this 
environment, dense star-forming gas has a clumpy structure 
which is off-center from the DM cusp. After powerful Pop 
III feedback ejects gas from this region, the Pop III remnants 
become less bound, further increasing the apocenter of their 
orbit. On the other hand, Pop II stars are rapidly produced 
only once the galaxy becomes relatively massive and concen- 
trated, and the DM center is more closely associated with the 
location of the densest gas. Pop II thermal feedback and ra- 
diative feedback is not strong enough to shut off further Pop 
II star formation near the galactic center, ultimately leading 
Pop II stars to have a very high spatial concentration. 

3.3. Epoch of Equivalence 

The dura tion of the Pop III epo ch has been previously stud- 
ied (e.g. iTrenti & Sti avelli 2009b by comparing the SFRs of 
the two populations, and finding the epoch where the Pop III 
SFR drops off dramatically compared to Pop II. Our Figure 
[2] shows that by a = 0.07 (z = 13), Pop II star formation is 
more common in all boxes except the void represented in Box 
UnderDense-. However this approach may not properly ac- 
count for the impact Pop III stars have on the universe at early 
times, as Pop III stars are able to generate significantly more 
thermal, ionizing, and metal feedb ack per baryon compare d 
to their Pop II counterparts (e.g. Tum linson & ShulH l2000). 
To more precisely quantify the impact of Pop III stars, we 
compute the relative "budget" of thermal energy, metals, and 
ionizing photons contributed in each galaxy by the two stellar 
populations, integrated over the lifetimes of all stellar parti- 
cles. For the purpose of this analysis, we define the "epoch 
of equivalence" as the time when Pop II stars have generated 
just as much ionizing radiation as their Pop III counterparts. 
In turn, we define the "duration of the Pop III phase" as the 
length of time from the formation of the first Pop III star to 
the epoch of equivalence. 

A global picture of the feedback budget is shown in Figure 
|4] Here, we calculate at each epoch the relative contribution 
of Pop III stars to the budget of ionizing photons, thermal 
energy from SNe, and metals injected into the ISM in run 
OverDense_nHle4_fid. Every budget can be examined indi- 
vidually. Pop II stars have produced as many ionizing photons 
as Pop III stars by a = 0.0633 (z = 14.8). The total thermal en- 
ergy released by SNe in Pop II stars surpasses that produced 
by Pop III stars at a = 0.0642 (z = 14.6). The total mass of 
metals produced in Pop II stars surpasses that of Pop III at 
a = 0.0689 (z = 13.5). Though we have chosen to define the 
epoch of equivalence in terms of ionizing radiation, any one 
of these feedback quantities leads us to the conclusion that the 
duration of the Pop III phase is short. 

To understand the transition to normal star formation in 
individual galaxies, we also budget the feedback from stars 
within a given galaxy. All numbers quoted in this section re- 
flect the integrated total feedback contributions over the life- 
time of stellar particles up to the epoch considered. As a 
case study, we take the third galaxy to form a Pop III star in 
run OverDense_nHle4_fid. The first star forms at a = 0.055 
(z = 17.2), when the halo mass and gas mass of this galaxy are 
9.8 x 10 6 M Q and 1.2 x 10 6 M Q , respectively. At a = 0.06 
(z = 15.7, 23 Myr later) this galaxy has formed a total of 
118M Q of Pop II stars. The first and only Pop III star that 
formed in the galaxy was a 17OM PISN progenitor. Al- 
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Figure 4. The fractional contribution of Pop III stars to the integrated feed- 
back budget of the universe in run OverDense_nHle4_fid vs scale factor. 
This fraction is calculated separately for metals (blue), ionizing radiation 
(black), and thermal radiation due to SNe (red). 



though the total mass of each stellar population is compara- 
ble, the one Pop III star still has had dominant influence on the 
galaxy by way of its stronger feedback. The Pop III star con- 
tributed 98.4% of the metals, 87.6% of the ionizing photons, 
and 93.4% of the thermal energy. This galaxy is clearly in the 
phase of its evolution where the effects of Pop III star forma- 
tion are most likely to be seen. However, we again note that 
despite this, the galaxy has not been significantly perturbed 
by the Pop III star's injection of energy. Only 1.7% of its hy- 
drogen is ionized, and all of the metals are confined in the gas 
well within the virial radius. 

By a = 0.0675 (80 Myr after the first star formed), one 
more Pop III star has formed, giving 270 M Q of Pop III mass 
formed in the galaxy, but now the total Pop II mass has 
jumped to 1262M . This second Pop III star has a metal- 
licity of log 10 Z/Z Q = -4.7, suggesting it formed in gas that 
was significantly enriched, but below the critical metallicity 
of log 10 Z/Z Q = -3.5. Therefore, whether the star should be 
Pop III or Pop II is rather sensitive to the particular setup of 
the model and stochastic effects. Pop III stars still contribute 
90.7% of the metals, 60.5% of the ionizing photons and 71.3% 
of the thermal energy. We see that the energy contributions 
from Pop II stars are rapidly becoming roughly equivalent to 
that of Pop III stars. Interestingly, the metal contribution of 
Pop III stars remains dominant here, confirming the trend of 
Figure 2] that the large yield of metals produced by PISNe is 
the longest lasting contribution of Pop III star formation. 

By a = 0.075 (z = 12.3, 134 Myr after the first star formed), 
the epoch of Pop III stars is clearly over in this galaxy despite 
the fact it has merged with another galaxy, which also hosted 
a recent PISN progenitor Pop III star. The combined initial 
mass of Pop III stars is now 440 M Q while Pop II stars con- 
tribute 48000 M . The Pop III stars still contribute 26.3% of 
the metal budget, but a mere 4.8% of the ionizing photon bud- 
get and 8.0% of the thermal energy budget. By a = 0.08, the 
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Figure 5. The duration of the epoch of Pop III dominance vs the halo 
mass of host galaxies at a = 0.1 (z = 9). Plotted are galaxies in ran 
OverDense_nHle4_fid (black), run UnderDense+_nHle4_fid (red), and run 
UnderDense-_nHle4_fid (blue). More massive galaxies have shorter Pop III 
epochs. 
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Figure 6. The epoch at which the contributions of Pop III and Pop II 
stars to the cumulative ionizing photon budget are equal is vs. the halo 
mass of host galaxies at a = 0.1 (z = 9). Plotted are galaxies in run 
OverDense_nHle4_fid (black), ran UnderDense+_nHle4_fid (red), and run 
UnderDense-_nHle4_fid (blue). The epoch of equivalence is earlier in more 
massive galaxies. 



Pop III metal contribution has finally dwindled away to only 
constitute 8.4%. 

The epoch when Pop III stars have affected the evolution of 
the galaxy is thereby confined to between the time when the 
first star formed (around a « 0.055) and the time when the 
budget of energy and metals of Pop II stars begins to dom- 
inate over Pop III stars (a « 0.07). The window of Pop III 
dominance only corresponds to da m 0.015, or 100 Myr. 

Figure|5]shows the duration of Pop III stars in this galaxy as 
well as in others, plotted against halo mass at a = 0.1 (z = 9). 
We note that the length of the Pop III epoch is generally be- 
tween 20-200 Myr, confirming that the galaxy discussed in 
the above paragraphs is not unique. The figure also reveals a 
trend where the most massive galaxies exit the Pop III stage 
earlier. These galaxies form in the most biased regions, where 
the infall of fresh gas overcomes the negative feedback gener- 
ated by the thermal and ionizing radiation of Pop III stars, and 
Pop II star formation erupts in the newly metal-enriched gas. 
Galaxies with higher mass at the time of Pop III star formation 
can also radiate the PISN energy away efficiently. 

Figure [6] shows that the epoch of equivalence for individ- 
ual galaxies is also strongly correlated with the mass of the 
halo at a = 0.1 (z = 9). This trend is perhaps not surprising, 
as the sites of formation for the first stars have to be the most 
biased regions of the simulation box, which would also pro- 
duce the most massive galaxies. It is impressive, however, 
that the trend is so definitive that it suggests we can predict 
the epoch at which Pop III stars became subdued based on the 
halo mass. 

We have also calculated Pop III epoch durations for galax- 
ies by using the metal and thermal energy budgets instead of 
the ionizing photon budget. The thermal energy dominance 
of Pop III stars tends to last on average about -10 Myr longer 
than the ionizing photon epoch, with little variance, while Pop 



III metals remain dominant for -20-50 Myr longer. These dif- 
ferences are significant, and can be understood by the relative 
efficiency of metal production per 170M Q Pop III star com- 
pared to Pop II counterparts. While our Pop II stellar parti- 
cles only generate an amount of metals equivalent to 1.1% of 
their initial mass, 170 M© Pop III stars release almost 47% of 
their mass in metals following a PISNe, meaning that Pop III 
metal feedback is -40 times more effective. In comparison, 
the number of ionizing photons per stellar baryon per lifetime 
is 6,600 and 34,500 for Pop II stars and 17OM Pop III stars 
respectively, translating to only a factor of -5 difference in 
efficiency. The variation between individual galaxies can be 
caused by the number of non-enriching 1OOM Pop III stars 
that form there, and by the rate at which Pop II stars form, 
as their feedback is released gradually over a longer interval 
of time. Despite this variation associated with using different 
feedback budget quantities, our conclusion on the brevity of 
the Pop III epoch within each galaxy remains robust. 

3.4. Extreme Pop III Feedback 

The duration of the phase of Pop III dominance may de- 
pend on the details of the implementation of stellar feed- 
back. Here, we consider two additional runs with arti- 
ficially enhanced Pop III feedback. By a = 0.0625 (z = 
15), run OverDense_nHle4_fid has 34 Pop III stars, while 
run OverDense_ExtremeSN has 28 stars and the Over- 
Dense_ExtremeRad run has 19 stars. At least part of 
this variation can be explained by the stochastic effect that 
the OverDense_ExtremeRad run happened to form 17OM 
PISNe progenitors more frequently than the other runs, hence 
quenching further Pop III star formation in their hosts. How- 
ever, other evidence suggests that the extreme radiation 
does indeed quench star formation very efficiently. This 
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can be discerned by examining Pop II star formation, as 
20000 M Q formed in run OverDense_nHle4_fid, 7000 M Q 
formed in run OverDense_ExtremeSN, and 4000 M© in run 
OverDense_ExtremeRad at a = 0.0625. While we showed 
in Paper I that the effects of baryon expulsion in individ- 
ual galaxies are temporary in both extreme feedback runs 
as well as the fiducial run, the cumulative effect of many 
more extreme feedback events appears to have a global ef- 
fect on the suppression of star formation at early times. More 
specifically, the relatively high IGM temperature in run Over- 
Dense_ExtremeRad (discussed in Paper I) means that newly 
accreted gas in this run will take longer to cool and condense 
to the densities required for star formation. It is apparent 
that in our code, radiative feedback is more effective than SN 
thermal feedback as a means of regulating star formation and 
modifying the structure of the ISM and IGM. 

In both extreme feedback runs, the epoch of equivalence 
for the relevant type of feedback is pushed back relative to the 
fiducial run, however it is always reached before z = 13. The 
epoch of equivalence for ionizing radiation increased by 5a = 
0.0073 relative to fiducial in run OverDense_ExtremeRad, 
while the epoch of equivalence for SN thermal energy in- 
creased by 5a = 0.0068 in run OverDense_ExtremeSN. Even 
drastic modifications to Pop III feedback output causes only 
marginal shifts in the duration of their epoch of dominance. 

3.5. Low Mass Pop III IMF 

Using a standard iMiller & Scalol (fl97l IMF for Pop III 
stars causes the duration of the Pop III epoch to be some- 
what shorter than the top-heavy fiducial case. The epoch 
of equivalence for ionizing radiation is reached in run Over- 
Dense_LowMass at a = 0.0622, a mere 90 Myr after the first 
star forms. 

However, the more important distinction in run Over- 
Dense_LowMass is the relative lack of suppression of star for- 
mation within individual galaxies. By a = 0.065 (z= 14.4), the 
total stellar mass formed in run OverDense_LowMass is a fac- 
tor of ~3 larger than in the fiducial run OverDense_nHle4_fid. 
The majority of stellar mass in both runs is contained in the 
three most massive star-forming galaxies, each of which has 
more stellar mass in the run OverDense_LowMass. In this 
run, nearly all of the stellar particles that formed in these three 
galaxies are within 10 parsecs of their galaxy's center, enrich- 
ing the dense gas slowly and steadily to perpetuate further 
star formation. A similar Pop II core is seen in the fiducial 
run, however the Pop III stellar remnants are often between 
10 and 70 pc from the galactic center. 

4. DISCUSSION AND CONCLUSIONS 

We find that Pop III star formation can be expected to ex- 
tend until z ~ 10. However, Pop III stars are not the dom- 
inant sources of any form of global feedback past a ss 0.07 
(z ~ 13) for our adopted cosmology, and only dominate indi- 
vidual galaxies for 20 to 200 Myr. 

The ionizing radiation contribution of Pop III stars is insuf- 
ficient to constitute a significant fraction of the photons re- 
quired to reionize the universe. Though they can temporarily 
evacuate gas, creating the opportunity for ionizing flux to es- 
cape the galaxy, it is unlikely that many photons were able 
to escape during these windows. We conclude this because 
in Paper I, we showed that accretion from filaments was able 
to bring back most of the expelled gas on roughly the same 
timescale (-150 Myr) that star formation in the galactic center 



was suppressed. This means that by the time star formation re- 
sumes at the center of the galaxy, the entire galaxy will again 
be surrounded by neutral gas, making it difficult for photons 
generated in the galaxy to escape to the IGM. 

We also showed in Paper I that most of the metals gener- 
ated in PISN would eventually fall back into the galaxy which 
hosted the PISN rather than remaining in the IGM. This fact, 
combined with the relative insignificance of the Pop III metal 
budget leads us to conclude that Pop III stars did not play 
much of a role in enriching the IGM, and instead that enrich- 
ment happened by normal stellar populations during and af- 
ter the epoch of reionization. Though Pop III star formation 
is self-limiting in individual galaxies, the absence of univer- 
sal enrichment implies that Pop III stars could form in un- 
derdense regions long after the universe has primarily tran- 
sitioned to normal star formation. It is in the least massive 
galaxies that formed in underdense regions where we might 
expect to find the most recent signatures of Pop III stars, sug- 
gesting once more that ultra-faint dwarf spheroidal galaxies 
may be the best place to look. As the chemical signature of 
PISNe metals is expected to be different from that of Type II 
SNe, and the metal budget of Pop III stars is the last feed- 
back tracer to be overtaken by Pop II, it would be prudent to 
continue searching for st ars that bear strange chemical abun- 
dances in these galaxies (Frebel 201 1). 

According to mock observational analysis for JWST 
by iPawlik et all d2012l) , galaxies with SFR in excess of 
0.1M Q yr" 1 may be detectable at z > 10. Only two galaxies 
in our simulation box meet this criterion in the final snapshot 
at z = 9, and both are evolved far beyond the point where Pop 
III stars made a significant contribution to their feedback bud- 
gets. A more promising way to detect galaxies that are still 
dominated by Pop III stars is by looking for PISN events di- 
rectly. This idea has been explore d and show n to be feasible 
for bo th LSST (iTrenti et al]|2009h and JWST dHummel et alJ 
|20T1 IWhalen et al.ll2012h A gain, considering our findings 
that underdense regions of the universe may stay pristine un- 
til late times, we conclude that the best approach for observ- 
ing campaigns is to focus on areas where the spatial clustering 
signal is low, such as cosmic filaments that connect two bright 
galaxies, to maximize the number of bright, low-redshift de- 
tections. 

If Pop III star s instead formed with a normal IMF (e.g. 
iGreif et aT] 1201 ll) . they would have even less impact on the 
universe, as demonstrated by our low-mass IMF run. While 
our extreme feedback runs showed that the Pop II transi- 
tion could be delayed with stronger Pop III feedback param- 
eters, the transition inevitably happened before z = 13 and 
takes less than 250 Myr for even the least massive galax- 
ies. In both cases, our modifications to the feedback out- 
put were simplistic ways to test the range of potential im- 
pact for Pop III stars without detailed modeling of additional 
physical processes. However, it is becoming clear that miss- 
ing physics in galaxy formation simulations can have a wide 
variety of effects that are not accounted by our limited ap- 
proach. For example, momentum-driven winds from radi- 
ation pressure can significantly stir the ISM increasing its 
susceptibility to o ther forms of feedback (I Wise et al.l 2012b; 
Agert zet alj2012l) . The presence or absence of dust can affect 
the shielding capabilities of molecular clouds from UV ra- 
diation, thereby modifying fragmentatio n and potentially the 
IMF of stars (lAvkut alp & Sp aansll201 ll) . Dark matter anni - 
hilation dSmith et al.ll2012l) . cosmic rays dUhlig et all 120 121) . 
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magnetic fields dTurk et aT1l2012l). and X-ray feedbac k from 
the first black holes (lHaimanll201 U Ueon et all 120121) could 
also be potentially important processes. We plan to explore 
them further in the next generation of simulations. 
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